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We investigate the low temperature behavior of a heavy 
particle (mass M) in a fermionic bath. An effective action is 
considered which exactly implements the orthogonality catas- 
trophe. It is equivalent to a model of local bosons coupling 
via a long-range interaction to the particle. The saddle point 
solution for the heavy particle propagator yields a finite quasi- 
particle weight in two and higher dimensions. In one dimen- 
sion the propagator decays algebraically for long times. The 
saddle point approximation is valid for large M. 
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The low temperature (T— >0) behavior of a heavy par- 
ticle strongly coupled to a fermionic bath has been a 
fascinating puzzle for more than two decades and yet has 
not been solved. Specifically, the question has been raised 
if, for T ~ >0, the particle has a plane wave character, or 
if the particle becomes a totally incoherent excitation so 
that the quasi-particle (QP) picture fails. The low tem- 
perature mobility of the heavy particle and the existence 
of a finite Drude weight are related problems. 

A variety of real physical systems are assumed to be 
'heavy-particle-fcrmionic-bath' systems: e.g. ions in nor- 
mal liquid 3 He diffusion of myons and hydrogen 
in metals ||, photocmission spectra of metals [||, and 
valence band holes in n-type doped semiconductors ||. 
Apart from these applications, a genuine theoretical in- 
terest in the solution of such a principal problem exists: 
it is closely related to the Kondo problem and the dis- 
sipative 2-level system |?]]. Further, it may give insight 
into the possible breakdown of Fermi liquid theory for 
strongly correlated electron systems: a single overturned 
spin in an otherwise spin-polarized Fermi gas with local 
interactions (Hubbard model) may or may not display a 
QP behavior |8p|l2|]. The existence of a QP-peak in the 
spectral function of the overturned spin will depend on 
dimension — and possibly on its effective mass and inter- 
action with the Fermi bath. Of course, the (non-) exis- 
tence of a QP-peak in the single-spin flip case does not 
predecide about the (non-) Fermi liquid behavior for the 
paramagnetic case. But it may be assumed that a Fermi 
sea of flipped spins will rather stabilize the Fermi liq- 
uid behavior because of the restricted phase space avail- 
able for scattering processes. Also, the overturned spin 
is not necessarily a heavy particle. However the case of 
a heavy particle should be investigated first because it 
is the most probable to reveal non-QP behavior, as the 
infinitely heavy particle does. 

In this paper we will focus on the question if the parti- 



cle shows coherent behavior for T— >0, i.e. if its spectral 
function displays a 5-function peak for zero momentum. 
The physics of this low-energy sector of the system is de- 
termined by two competing phenomena: orthogonality 
catastrophe (OC) [13[ versus recoil |JT§]. On the one 
hand, we expect an incoherent behavior because of An- 
derson's OC valid for an infinitely heavy particle: a sud- 
den displacement of the particle will produce an infinite 
number of low-energy particle- hole (ph) excitations. Sim- 
ilarly, a localized valence band hole generated through 
optical absorption at time will produce a shakeup of 
the conduction band in the long time limit. This causes 
a vanishing overlap of the original state with the final 
state in the thermodynamic limit. The final state is char- 
acterized by finite phase shifts of all the bath particles, 
necessary to screen the localized hole. On the other hand, 
momentum conservation for a finite-mass particle implies 
a finite recoil energy up to .©recoil — (2&/) 2 /2M of the 
heavy particle with mass M . This finite E leco {\ restricts 
the available phase space for scattering processes on the 
Fermi surface simply by energy conservation. Accord- 
ingly, a perturbative calculation in the interaction of the 
heavy particle with the fermionic bath shows that the 
number of low lying ph-excitations is reduced in space 
dimensions d > 1 so substantially that the QP-peak is 
restored for finite mass |J. However, it has been argued 
that conclusive evidence for a finite QP-peak in 2d has 
to be obtained from a non-perturbative approach — even 
if the interaction is small. 

To set up the problem we consider the following Hamil- 
tonian: 



H = H n + H] 



bath 



H, 



(1) 



Here, the kinetic energy of the heavy particle with mo- 
mentum P is 

H = P 2 /(2M) 

For some specific lattice cases, it would be more appro- 
priate to assume a band dispersion. We will comment on 
such a case below. Actually, it is important to include 
a periodic potential for the discussion of self-trapping 
|l5| Il7| . But it has been asserted that, even for charge 
Z > 2, localization is unlikely to occur Jl8| |. 

The bath contribution is just a Fermi (lattice) gas with 
dispersion e& 
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Spin may be included, and the results have to be modified 
in a trivial way by a spin degeneracy factor. The inter- 
action is assumed to be local, i.e. it acts at the position 
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R of the heavy particle: 

Higher scattering channels beyond s-wave should be con- 
sidered for an 'extended particle' but here we will inves- 
tigate only s-wave scattering to keep the mathematics 
transparent. 

We will present a path integral approach which cor- 
rectly reproduces the known limits (in any dimension): 
the second order perturbation theory in U |l9|], the high 
temperature case (-E roC oii < T <C ef), and the localized 
case (M = oo) which straightforwardly results from the 
high temperature case (with -Brecon = 0). As regards the 
coherence of the heavy particle state, we find: algebraic 
decay of the heavy particle propagator in lei, but finite 
QP-weight for higher dimensions. Hence, the finite QP- 
weight in 2c? is the major result of our considerations. 

The partition function for (0) is Z = Z Q fVRe~ lS W . 
Z is the partition function for the fermionic bath. S — 
S + Su splits into the action for the free particle S — 
\ Jq drM R(t) 2 and the coupling to the bath 

S v = trln[l + {/ 5fc (r-r')e j(fc - fc ' ) - R(T ' ) ] (2) 

which results from the integration over the fermionic 
fields. <?i,(t) is a free-fermion propagator. The argument 
of the logarithm is a matrix indexed by r, r', fc and k'. 
Its power series will generate terms involving any num- 
ber of time integrals. In the following paragraphs we will 
argue that (||) may be cast into a 'singly retarded' form 
which requires two time integrations, only. 

A path integral which builds on all the (infrared) singu- 
lar contributions from paths connecting two heavy parti- 
cle positions was, to our knowledge, first proposed by Sols 
and Guinea |L7| . Although we have to refer the reader to 
Prokof 'ev's papers |^0| for the detailed discussion of the 
set-up of the action, we want to highlight the connection 
to the orthogonality catastrophe (OC). The OC not only 
arises for the overlap of a state with undisturbed Fermi 
sea and a state with all particles scattering from a lo- 
calized impurity. But also the overlap of the electronic 
ground state 0^ with impurity at position R and a state 
(j} R > vanishes in the thermodynamic limit N e — > oo. N e 
is the number of electrons (bath particles). Yamada and 
Yoshida et al. managed to crack this integral: 

(^ R )^e^ R '- R ^^ (3) 

where 



F(\R\) = s arcsrnWl - x 2 sin 5) (4) 

for pure s-wave scattering with phase shift <5, and 
x = (e R )\k\ = kj, which is cos(kfR), J a (kfR) and 
sm(kfR)/(kfR) for 6?= 1, 2 and 3, respectively. F(R) 
is also known for extended scatterers with phase shifts in 
higher order angular momentum channels [p2|. Now we 



consider the general form of a singly retarded (self-) in- 
teraction. Due to translational invariance it is certainly a 
functional of R(t')—R(t), and also of t'—t. We are only 
interested in the low frequency behavior which restricts 
the scattering processes to the vicinity of the Fermi sur- 
face where energy and momentum are independent vari- 
ables. Accordingly, the retarded interaction separates 
into a space and time dependent part 0,|o) : 

S cS = S - X - JJ^dTdT'T{\R'-R\) P {t'-t) (5) 
p is written as 

which reflects the fact that the number of excitations 
close to the Fermi surface is proportional to uj, and cor- 
rections in to 2 are irrelevant for the investigated infrared 
divergence. Evaluation of the overlap (^) and comparison 
with (§ fixes T: 

T{R) = (-) 2 - F(\R\) 

which is displayed in the figure. It includes all paths 
responsible for the OC Eq. (|) @. 

The effective action (||) may be verified for high tem- 
peratures or for small interaction U by expanding (|^) 
either in R—R' J2l| or in U. The second order in either 
expansion is consistent with this singly retarded action 
S e g which, moreover, covers contributions from any or- 
der. 

There have been attempts in the past to extract infor- 
mation about the low-T behavior from this action [ p"7| ]. 
Either it is assumed that the quadratic terms in R already 
determine qualitatively the low-T behavior |pp| . Then 
the system is actually equivalent to a harmonic oscillator 
bath with spectral density ~ u> coupled to the heavy par- 
ticle (Caldeira-Leggett model M). Such a mapping onto 
a bosonic bath was shown to be correct for the X-ray edge 
problem with a localized particle |24| . Another approach 
mimics the shape of F{R) for d > 2 by assuming a step- 
like functional form with the step at R equal the lattice 
constant [ p5| . The latter approach is equivalent to a pro- 
cedure by Hedegard and Caldeira [^6| which introduces 
independent oscillator baths at each lattice site. They 
further assumed that if the particle comes back to a lat- 
tice site twice or more times it always encounters a "new" 
bosonic bath in its thermal equilibrium, i.e. the history 
of excitations at a site is not recorded ('non-interacting 
blip approximation'). Again, as for the harmonic oscilla- 
tor bath, no QP-peak shows up for T — > 0. 

We will demonstrate that the details of the long- 
distance behavior of T{R), i.e. the Friedel-oscillations, 
actually control the spectral function of the heavy parti- 
cle. Since a direct integration over the field R is hope- 
lessly complicated we propose the following detour: to 
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disentangle the untractable retarded part of the action 
we introduce new bosonic fields 6 7 fc(r). Then the prob- 
lem is equivalent to the following polaron-type model, as 
can be checked by integrating out the bosons: 

H b = H + J2 E 7 b\ h hy h + ]T(^ b\ k J hR + c.c.) (6) 

yk yk 

with a spectral density X) 7 ~ ^7) = w0(cj) for 
the local bosons, and a long-range interaction l^) 2 = 
J2r F(R)e lkR . The particle coordinate is eliminated by 
a transformation [pTj Tr — exp[i(^ 7fc k b* lk b lk ) R ] on 
the local coordinate system of the heavy particle. The 
transformed Hamiltonian reads 

H Po = H khl {P ) +^2 E f h \k h ik + J2( y k b \k + c - c -) 
with 

H kin (Po) = (Po -Y / kb \k b ik) 2 /(2M) 
7 fe 

where the c-number P a is the total momentum of the 
"polaron" . With this transformation, V k is just a dis- 
placement of the harmonic oscillators whereas the ki- 
netic energy of the heavy particle introduces a compli- 
cated current-current interaction, the origin of which is 
the finite recoil. 

We are interested in the real-time heavy particle prop- 
agator (at T = 0): 

iG Po {t) = (0|(P o | e -^ 4 |P o )|0) ~[v{b*,b]e- lSp ° lb * Mt] 



with 



|0) is the Fermi sea which then translates to the vacuum 
for the bosons of (^|), and Sp o is the action corresponding 
to H Po : 

S Po = S {P ) + V f dt(V k b; k + v: 6 7fc ) + Si n t 
_.«. Jo 



with 



/"* k 2 kP 

s o(Po) = E / dt Kk(-^t + E i + Yn-^r) b ^ 

yk J ° 



and 



S'int. — 



1 f* 

— J dtkb* k b jk -k'b*, k ,b jlk , 

yk.y'k' 



A saddle point evaluation of Gp o= o{t) yields b° k {t) = 

-iV k J*dt*e i( ^ E -' +k2 / 2M ^- t) and the action at the sad- 
dle point is: 

S° P=0 = -t + i f duJ ^l(e-^ - 1) 

J U) J u z 



Qd{") =E lk \Vk\ 2 6(Lo-E, ( ~£j) 



= fdn k j^dkk^\v k \ 2 ^-^) 

First, we consider the case of infinite mass in any di- 
mension: we find correctly £>^°(<^) = u>F(R = 0) — uj(-) 2 , 
i.e. the exponent for the algebraic long-time decay of the 
heavy particle propagator, G(t) ~ 9{t) \t\~ a , is a = (^) 2 



In d = 1, 3~{R) is a periodic function. Therefore its 
Fourier spectrum is discrete and only the k = compo- 
nent contributes to the infrared behavior so that 



q 1 (cj) = uj d(k f R)T{R) for uj -> 0, M ^ 00 
Jo 

The OC-exponent for the algebraic long-time decay is 
a = Cd{kfR)T{R), i.e. a = f(5)-(^) 2 with /(£-►£) = 
2/3 j2| and f(6-*Q) = 1/2. It is well understood that a 
for M^oo is different from a for M = oo [|ll[: the finite 
recoil for M^oo excludes scattering processes across the 
Fermi sea for u> < E Teco n, whereas these processes cost no 
energy for M = oo. 

Half-filling in d = 1 is special for a nearest neighbor 
hopping motion of the heavy particle. The oscillations 
of J-(R) are commensurate with the lattice periodicity 
which implies that a hopping heavy particle will probe 
!F{R) only at its maximum and minimum, whereas a 
'free' heavy particle averages over all values of T{R). 
Therefore we expect for this special commensurate case 
a = i((5/7r) 2 , even for S — > This compares well with 
numerical findings by H. Castella |Tj|j . 

In d = 2 



2 5 tan S 



,3/2 



y / fc 2 /2Af 



for 



0, M/oo 



i.e. the power in uj is larger than 1 so that the QP-peak is 
finite and the particle displays coherent behavior. The in- 
coherent part of the spectral function diverges for oj — > 
as 1/s/u, as already known from second order perturba- 
tion theory ||. The actual value of the QP-weight and 
the prefactor of the incoherent part are determined by 
the details of the high frequency behavior. 

Only if the Fermi energy is located at a van Hove sin- 
gularity 02^) becomes linear in w, as can be checked in 
the 2 nd -order perturbation theory. The QP-peak in 2d 
and its suppression at a van Hove singularity was also 
found numerically by S. Sorella |l2| using a variational 
ansatz p0| . 

For d = 3 



IS tan S 
Q3{u)= 7 — 

4 7T Z 



k 2 f /2M 



for 



0, M/oo 



which yields a QP-peak as in d - 
incoherent background. 



2, but on top of a finite 
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We further want to mention that the step function ap- 
proximation to T{R) will produce a finite QP-peak in 
any dimension: 

The incorrect frequency dependence of g s t cp illustrates 
that the details of the long distance behavior of T{R) are 
indeed essential. This result has to be contrasted with 
that of the non-interacting blip approximation which is 
derived by skipping the term k 2 /2M in S : the QP- 
weight is zero in any dimension due to the missing recoil. 

To estimate corrections to the saddle point result we 
calculated the fluctuations at the stationary point and 
the leading order shift of the saddle point by fluctuations. 
All corrections are of 0(1/M). Therefore, the question 
may be raised if the result for the Id spectral function 
is exact for large M and small w (lo <C -E^ccoU *C 
We expect the result for d — 1 to be quantitatively cor- 
rect if the limit of heavy mass can be implemented by a 
singly retarded action. The spectral function in d > 1 can 
be only qualitatively valid since the (non-zero) value of 
the QP-weight and the prefactor of the incoherent back- 
ground depend on the details of the high energy behavior. 

We acknowledge many helpful discussions with 
N. Andrei, H. Castella, F. Mila, E. Muller-Hartmann, 
N.V. Prokof'ev, A.E. Ruckenstein, and P. Wolflc. This 
work was supported by the DFG (T.K.). 



[14] E. Muller-Hartmann, T.V. Ramakrishnan, G. Toulouse, 

Phys. Rev. B 3, 1102 (1971). 
[15] A. Schmid, Phys. Rev. Lett. 51, 1506 (1983). 
[16] M.P.A. Fisher, W. Zwerger Phys. Rev. B 32, 6190 (1985). 
[17] F. Sols, F. Guinea, Phys. Rev. B 36, 7775 (1987). 
[18] D.S. Fisher, A.L. Moustakas, to appear in Phys. Rev. B 

51, (March 15, 1995). 
[19] as given in a cumulant expansion. 

[20] N.V. Prokof'ev, J. Mosc. Phys. Soc. 2, 157 (1992); Int. 

J. Mod. Phys. B 7, 3327 (1993). 
[21] K. Yamada, K. Yoshida, Prog. Theor. Phys. 68, 

1504 (1982); K. Yamada, A. Sakurai, S. Miyazima, 

H.S. Hwang, Prog. Theor. Phys. 75, 1030 (1986), and 

references therein. 
[22] K. Vladar, Prog. Theor. Phys. 90, 43 (1993). 
[23] R(t) — const is the high-T saddle point. 
[24] K.D. Schotte, U. Schotte, Phys. Rev. 182, 479 (1969). 
[25] A. Rosch, diploma thesis (Univ. Karlsruhe, Jan. 1994). 
[26] P. Hedegard, A.O. Caldeira, Physica Scripta 35, 609 

(1987); P. Hedegard, Phys. Rev. B 35, 533 (1987). 
[27] T.D. Lee, F.E. Low, D. Pines, Phys. Rev. 90, 297 (1953). 
[28] P. Nozieres, C.T. de Dominicis, Phys. Rev. 178, 1097 

(1969). 

[29] S — > 7r/2 corresponds to U — > oo. 

[30] D.M. Edwards, Prog. Theor. Phys. Suppl. 101, 452 
(1990). 



FIG. 1. T{R) for d = 1, 2, and 3 with S close to tt/2. 



[1] B.D. Josephson, J. Lekner, Phys. Rev. Lett. 23, 111 
(1969). 

[2] P. Wolfle, in: Lecture Notes in Physics, vol. 115, ed. by 
A. Pekalski, J. Przystawa (Springer, Berlin, 1980), and 
references therein. 

[3] J. Kondo, T. Soda, J. Low Temp. Phys. 50, 21 (1983). 

[4] Y. Kagan, J. Low Temp. Phys. 87, 525 (1992), and ref- 
erences therein. 

[5] S. Doniach, Phys. Rev. B 2, 3898 (1970). 

[6] M. Calleja, A.R.Goni, B.S.Dennis, 

J.S.Weiner, A. Pinczuk, S. Schmitt-Rink, L.N. Pfeiffer, 
K.W. West, J.F. Mueller, A.E. Ruckenstein, Solid State 
Commun. 79, 911 (1991). 

[7] A.J. Leggett, S. Chakravarty, AT. Dorsey, M.P.A. Fis- 
cher, A. Garg, W. Zwerger, Rev. Mod. Phys. 59, 1 (1987), 
and references therein. 

[8] T.B. McGuire, J. Math. Phys. 6, 432 (1965). 

[9] T. Kopp, A.E. Ruckenstein, S. Schmitt-Rink, Phys. Rev. 
B 42, 6850 (1990). 
[10] D.M. Frenkel, Phys. Rev. B 46, 15008 (1992). 
[11] H. Castella, X. Zotos, Phys. Rev. B 47, 16186 (1993); 

H. Castella, to be published. 
[12] S. Sorella, Phys. Rev. B 49, 12373 (1994). 
[13] P.W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 



4 



